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■ - - . ABSTRACT 

I We report on the spectral analysis results of the neutron star, atoll type, low mass X-ray Binary 

. 4U 1636-53 observed by INTEGRAL and BeppoSAX satellites. Spectral behavior in three different 

psj ' epochs corresponding to three different spectral states has been deeply investigated. Two data set 

spectra show a continuum well described by one or two soft blackbody plus a Comptonized components 
with changes in the Comptonizing electrons and black body temperature and the accretion rates, which 
are typical of the spectral transitions from high to low state. In one occasion INTEGRAL spectrum 
shows, for first time in this source, a hard tail dominating the emission above 30 keV. The total 
spectrum is fitted as the sum of a Comptonized component similar to soft state and a power-law 
component (F = 2.76), indicating the presence of a non thermal electron distribution of velocities. In 
this case, a comparison with hard tails detected in soft states from neutron stars systems and some 
^ . black hole binaries suggests that a similar mechanism could originate these components in both cases. 

fS| ' Subject headings: accretion, accretion disks - gamma rays: observations - radiation mechanisms: 

l/^ , non-thermal - stars: individual: 4U 1636-53 - stars: neutron - X-rays: binaries 

\o 
o 
\o 

. 4U 1636-53 is a neutron star low mass X-ray binary (LMXB) classified as a atoll source (Hasinger & van der Klis 
1989), with an orbital period of 3.8 hr derived from the optical variability of its companion V801 Arae (Pedersen, van 
Paradijs & Lewin 1984) and at distance of 3.7-6.5 kpc (Fujimoto et al. 1988, Smale & Lochner 1992, Augusteijn et 
al. 1998). 

O While the X-ray burst properties and timing signatures have been analyzed extensively (see Jonker et al. 2005, 
Belloni et al. 2005 and references therein) the spectral characteristics have been studied only at low energy with 
CZJ Einstein , EXOSAT , Temna and ASGA . In general, the spectrum was acceptably fitted by a Comptonization model 
plus a black body component. A coronal temperature of 2.3-2.6 keV, an optical depth of 13, and a soft black body 
^ , temperature of 0.55 keV were the best fit parameter for the Einstein data (Christian and Swank 1997). EXOSAT 
• ■ data reveal two different source intensity: for the maximum intensity the coronal temperature was ~ 3.5 keV and the 
' optical depth ^ 10.4, while for the minimum intensity the temperature was ^ 6.1 keV and the optical depth ~ 7.2. 
H I White et al. (1988) found the coronal temperature to be 1.8 keV and the optical depth to be 16.7. A model consisting 
. - - 1 of a black body, a multi-color disk and a broad Fe line, was used by Asai et al. (1998) to fit the ASGA data, confirming 
previous Temna results (Waki et al. 1984). We report here a broad band spectral analysis performed on data from 
BeppoSAX and INTEGRAL satellites, which allowed us to better constrain the spectral parameters and to detect the 
presence of a high energy tail dominating the spectrum above ~ 30 keV. A similar feature has been observed in other 
LMXBs namely GX 17+2 (Di Salvo et al. 2000), GX 349-1-2 (Di Salvo et al. 2001), Sco X-1 (D'Amico et al. 2001), 
4U 1608-522 (Zhang et al. 1996), XB 1254-690 (laria et al. 2001), Cir X-1 (laria et al. 2002) and 4U 0614+091 
(Piraino et al. 1999). 

2. OBSERVATIONS AND DATA ANALYSIS 
Table El summarizes the log of INTEGRAL and BeppoSAX observations of 4U 1636-53 . 

BeppoSAX observed the source on three occasions: February and March 1998 and February 2000. LEGS, MEGS and 
PDS event files and spectra , available from th e ASI Scientific Data Center, were generated by means of the Supervised 
Standard Science Analysis ijFiore et al. 19991) . Both LEGS and MEGS spectra were accumulated in circular regions of 
8' radius. The PDS spectra were obtained with the background rejection method based on fixed rise time thresholds. 
Publicly available matrices were used for all the instruments. The cross-calibration constant values were taken in 
agreement with the indications given in Fiore et al. 1999. Fits are performed in the following energy band: 0.5-3.5 
keV for LEGS, 1.5-10.0 keV for MEGS and 15-70 keV for PDS. 

The analyzed INTEGRAL (Winkler et al. 2003) data set consists of all observations in which 4U 1636-53 was within 
the high-energy detectors field of view. Observations are organized into uninterrupted 2000 s long science pointing, 
windows (sew): light curves and spectra are extracted for each individual sew. Wideband spectra (from 5 to 150 
keV) of the source are obtained using data from the two high-energy instruments JEM-X (jLund et al. 200,11 and IBIS 
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TABLE 1 

BeppoSAX and integral Observations 



BcppoSAX Journal 
Start Date Exposure time Count s 









LEGS 


MEGS 


PDS 


LEGS 


MEGS 


PDS 








ksec 


kscc 


kscc 


[0.4-3 kcV] 


[1.5-10 kcV] 


[20-60 kcV] 




epoch (a) 


1998-02-24 


13 


39 


17 


14.66 ± 0.03 


39.37 ± 0.02" 


0.48 ± 0.04 




epoch (b ) 


1998-03-01 


6 


14 


7 


14.95 ± 0.05 


25.42 ± 0.04'' 


0.80 ± 0.06 


1»' 


epoch (c) 


2000-02-15 


12 


37 


19 


22.08 ± 0.04 


29.09 ± 0.03' 


0.73 ± 0.03 










INTEGRAL Journal 














Start Date 


Exposure 


time 


Count s ^ 












JEM-X 


IBIS 


JEM-X 


IBIS 












ksec 


kscc 


[5-15koV] 


[20-150 kcV] 








2""* epoch 


2003-03-04 


36 


594 


6.3 ± 0.2 


10.74 ± 0.08 








3^^ epoch 


2003-03-04 


16 


117 


3.70 ± 0.06 


3.22 ± 0.04 





" MEGS count rates refer to MEGS2 and MEGS3 units. 
MEGS count rates refer to only MEGS2 unit. 



ljUbertini et al. 2003j) . Data were processed using the Off-line Scientific Analysis (OSA version 5.1) software released 
by the INTEGRAL Scientific Data Centre. While IBIS provide a very large FOV 30°), JEM-X has a narrower 
FOV (~ 10°), thus providing only a partial overlap with the high-energy detectors. Data from the Fully Coded field 
of view only for both instrument have been used. The angular resolution of IBIS instrument is 12 arcmin, any source 
at a distance larger than the instrument angular resolution do not contribute at all to the observed source spectrum 
and flux (Ubertini et al. 2003) due to the coded mask intrinsic characteristics. 

3. SPECTRAL ANALYSIS RESULTS 

The whole data set was carefully fitted with several physical models, while trying to keep the number of free 
parameters as low as possible. 

Each time a new component was added to the model, a F-test was performed. We assumed that a F probability 
larger than 95% implies a significative improvement of the fit. The uncertainties are at 90% confidence level for one 
parameter of interest (Ax^ = 2.71). 

When spectra are from more than one detector, we allow the relative normalization to be free with respect to the 
MFCS and IBIS data, for BeppoSAX and INTEGRAL respectively. XSPEC v. 11.3.1. has been used. 
Spectral behavior has been studied separately in three epochs consisting of the following data: 

1** epoch: all three BeppoSAX observations available from February 1998 to February 2000. During these periods 
the source was always in a soft/high state. 

2"'' epoch: JEM-X and IBIS data available from 52644 MJD to 53644 MJD with count rate < 5 counts s^^ in the 
20-40 keV energy band. For the chosen period the source was in a hard/low state. 

epoch: JEM-X and IBIS data available from 52644 MJD to 53644 MJD with count rate > 5 counts s^^ in the 
20-40 keV energy band. This epoch does not corrspond to either the soft or the hard state and and here we call it 
peculiar state as will be explained in detail later on. 

The most simple model which provides a good fit to each BeppoSAX spectrum in the energy band 0.5-70 keV consist 
of a thermal Comptonized component modeled in XSPEC by gomptt (a spherical geometry was assumed) plus a soft 
component which we modeled with two temperature blackbody. The simplest model consisting of multicolor diskbb 
(Makishima et al. 1986) plus a gomptt component do not give a good fit, with a ~ 2.5 for each observations. 
The black body and thermal Comptonized component parameters were left free in order to determine the blackbody 
temperature T;,;,, the electron temperature Tg, the optical depth r and seed photon temperature Tq. Results from 
these fits are reported in Table [5] Figure ^ shows three BeppoSAX spectra and the residuals with respect to the 
corresponding best fits. The column density Nh towards the source was left free and its value measured by the LECS 
and MFCS instrument is always close to the galactic column density (A'h galactic = 3.58 x 10^^ cm^^, estimated 
from the 21 cm measurement of Dickey & Lockman 1990). The seed photon temperature is always To=1.3 ± 0.2 keV, 
and all parameters did not significantly vary from one observation to another. 

In order to achieve the highest signal to noise ratio we build the BeppoSAX average spectrum arranging all three 
observations. This procedure can, in principle, be risky since the source can change its spectrum from one observation 
to another. However in our case, the previous analysis showed no significant shape changes for the three observations. 
We then can take advantage of the high quality of the average spectrum up to about 70 keV. We fit the average 
BeppoSAX spectrum with the model used for the single observations. Nh has been fixed to the value for the galactic 
column density. Spectral fit results are given in Table |21 the average spectrum is shown in Fig.|3 

The source was in the soft/high state with an un-absorbed luminosity of Lq. 1-200 keV — 2.0 x 10'^'' erg s~^, assuming 
a distance of 5.9 kpc (Cornelisse et al. 2003). As in the case of the Einstein observation (Christian and Swank 1997), 
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TABLE 2 

Results of the fit of 4U 1636-53 spectra in the energy band 0.5 — 70 keV for three BeppoSAX observations. For each 

OBSERVATIONS WE REPORT THE BEST FIT MODEL. THE BLACK BODY EQUIVALENT RADIUS IS COMPUTED USING A DISTANCE OF 5.9 KPC 
(CORNELISSE FT AL. 2003). UNCERTAINTIES ARE AT THE 90% CONFIDENCE LEVEL FOR A SINGLE PARAMETER VARIATION. 



BeppoSAX spectra 

'model — bbody + bhody + comptt 

Nh TbBI TbB2 To T iiBBl RbS2 "COMPTT X^/d.O.f 

10^^ keV keV keV km km 

3.1^0 4 0.28 ± 0.03 0.61 ± 0.03 3. 4t°J sfl 57^1^ 2Atl 0.wto%l 190/124 

3.4+°;^ 0.26 ±0.04 0.59+°;°^ 2.9+° j btl Sltf^ 23l| 0.241° °^ 132/122 

3.2to 2 0.28 ± 0.03 0.60 ± 0.03 2.8^04 5±2 &9>t\% 21 ± 2 O.ieton^ 133/122 



o 
o 




o — 

in - 




^ , — , — ^ , , , , , — , — I 

1 10 
channel energy (keV) 



Fig. 1. — Three BeppoSAX spectra of 4U 1636-53 and the residuals with respect to the corresponding best fits, in the range 0.5-70 keV. 



thermal Comptonization of the optically thick plasma corona with a quite low electron temperature is dominating. 
The most simple model which provides a good fit to the INTEGRAL hard state consist of a thermal Comptonized 
component modeled in XSPEC by comptt (Titarchuk 1994) (a spherical geometry was assumed) plus a soft component 
which we modeled by a single temperature blackbody. Because of the very good low energy BeppoSAX coverage, we 
used the average value of input soft photon temperature and column density (Tq = l.SkeV and Nh = Nh Galactic) 
measured by the BeppoSAX observations. During this period, the source has a luminosity of io. 1-200 kcV — 1-4 x 10^'' 
erg s~^, lower than the one in the soft state, in agreement with the usual ranking of the luminosity in atoll sources 
(e.g., Hasinger & van der Klis 1989; van der Klis 2000; Gierlinski & Done 2002). The electron temperature is now 
substantially higher than in the soft state. To ~ 23 keV, and the Comptonization component extends well above 100 
keV. 

The same model was been applied to the 2'"* INTEGRAL data set but resulted in a poor fit with a x^/d-O.f = 97/57 
and clear residuals above 60 keV. Adding a power law improves the fit significantly {x^/d.o.f becomes 64/55). A 
simple model with the single temperature blackbody and power law do not fit our data {xlid 5). The disk 
component becomes negligible and it is not necessary to best fit the data. Spectral fit results are given in Table |21 and 
spectra are shown in Fig. |31 . In this peculiar state, the spectrum is also compatible with a power-law and multiple 
blackbody (Makishima et al. 1986) instead of Comptonization halo; nevertheless the temperature is extremely high 
and perhaps unphysical. Since the spectral parameters in the peculiar state are similar to those in the soft/high 
BeppoSAX state, a possible interpretation of this state is that 4U 1636-53 was in a similar soft/high state as during 
the BeppoSAX observations, but with a new overlapped component at high energies simultaneously present. The 
un-absorbed luminosity is io.i-20okcV — 7.4 x 10'^^ erg s~^. Comparison between the models in the three different 
spectral states are shown in Fig. 01 Spectral state transitions are evident, with the hard INTEGRAL state extended up 
well above 100 keV and, in the peculiar state, a steep power-law component detected with high statistical significance 
up to 100 keV. 



4. DfSCUSSfON AND CONCLUSfONS 

According to our present understanding, the black-body component in the soft states could originate at both the 
neutron star surface and the surface of an optically-thick accretion disk. In our observations the two black body 
components seems to originate from two different parts of the disk, corrisponding to two different temperature. The 
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TABLE 3 

Results of the fit of 4U 1636-53 spectra in the energy band 0.5 — 70 keV and 5 — 150 keV for BeppoSAX and INTEGRAL 

OBSERVATIONS RESPECTIVELY, FOR THREE DIFFERENT SPECTRAL STATES. FOR EACH STATE WE REPORT THE BEST FIT MODEL. THE BLACK 
BODY EQUIVALENT RADIUS IS COMPUTED USING A DISTANCE OF 5.9 KPC (CORNELISSE FT AL. 2003). UNCERTAINTIES ARE AT THE 90% 

CONFIDENCE LEVEL FOR A SINGLE PARAMETER VARIATION. 





BeppoSAX average high/soft state spectrum 






model — bbody + hi 


ody + comptt 






Tbbi TbB2 


To T„ T 




^BB2 


racOMPTT X^/d.o.f 


keV keV 


keV keV 


km 


km 




0.24 ± 0.01 0.58 ± 0.02 


1.3 ± 0.1 3.4 ± 0.3 3.8 ± 0.4 


79 ± 6 


25 ± 1 


0.19 ± 0.02 235/190 




INTEGRAL low/hard state spect 


rum 






model — hbody + comptt 






TbB2 


To Te T 




^^COMPTT 


xVd.o.f 


keV 


keV keV 


km 


10"^ 




1 9 + 0.2 
-^■^-0.3 


1.3 fixed 23^2 1-lto 3 


r + 3 
^^-2 


1 0+"-^ 


57/59 



INTEGR,AL peculiar state spectrum 

model — comptt + powerlaw 

To To T r nCOMPTT "powerlaw X^/d.O.f 

keV keV 10"^ ph keV'^ cm^^ 



1.3 fixed 5^3 3±2 2.6 ±0.1 2.8l!;-g ^-^-ot '''2/57 




channel energy (keV) 



Fig. 2. — The BeppoSAX average spectrum in the soft state (epochs 1), shown together with the total model and its components: the 
total, two blackbody and the Comptonization components are shown in green, red, magenta and blue, respectively. Residuals with respect 
to the corresponding best fits are also showed. 



high optical thickness of the plasma rules out the black body component as originates from the neutron star cap 
regions. Anyway, independently from the optical thickness of the plasma, it is very improbable to have emission from 
the polar caps of the neutron star for this source. Infact the magnetic field of the LMXBs is ~ 10* Gauss implying 
that the accretion onto the neutron star is not magnetically dominated and the matter accretes onto the whole NS 
surface. The Comptonization component may arise from a corona above the disk and/or between the disk and the 
neutron star surface. In the hard states, accretion probably assumes the form of a truncated outer accretion disk and 
a hot inner flow, joining the disk and the stellar surface as previously reported by Barret & Olive (2002) for LMXB 
4U 1705-44. The spectral transitions are generally, but not necessarily, coupled with changes in luminosity, indicating 
they are driven by variability of the accretion rate or change of the geometry of the system as for Black Hole hard/soft 
transition at constant luminosity (Belloni et al. 2005). For 4U 1636-53 the accretion rate is lower in the hard state 
than in the soft state: for an accretion efficiency of 77 = 0.2 (corresponding, e.g., to Af^g = 1.4Mq and i?NS = 10 
km ) and using our model luminosities, io. 1-200 kcVi we get Msoft — '7.3 x 1O~^M0 yr^^ and Mhard — 4.0 x lO^^M© 
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Fig. 3. — The spectra of epochs 2 (soft state, left) and 3 (peculiar state, right) observed by INTEGRAL , shown together with the 
total model and its components. Left: the total, the blackbody and the Comptonization components are shown in green, magenta and 
blue, respectively. Right: the total, the Comptonization and the power law component are shown in green, blue and orange respectively. 
Residuals with respect to the corresponding best fits are also showed. 




Energy (keV) 

Fig. 4. — Comparison of the models for the three observations. Epoch 1, BeppoSAX soft state: green; epoch 2, INTEGRAL hard state: 
red; epoch 3, peculiar INTEGRAL state: blue. 



yr ^. Transitions are apparently accompanied by changing in the geometry of the flow, and in the relative contribution 
of the blackbody-like and Comptonization components. In the peculiar state the accretion rate becames very high 
{Mpecuiiar — 2.1 X IO~^Mq yr~^), but this value can be influenced by the steep power law that dominates the energy 
spectrum at low energies. The emission is best described by Comptonization from a complex electron distribution due 
to a low temperature (^^ 3 — 4keV) thermal electron distribution together with non-thermal power-law electrons. This 
two-component electron distribution could arise from non-thermal electron acceleration regions powered by magnetic 
reconnections above a disc. Low-energy electrons cool preferentially by Coulomb collisions leading to a thermal 
distribution while the high-energy electrons cool by Compton scattering, preserving a non-thermal distribution (Coppi 
1999). Alternatively, the thermal and non-thermal electrons could be spatially distinct, e.g. magnetic reconnection 
above the disc can produce a non-thermal electron distribution, while overheating of the inner disc produces the 
thermal Comptonization (Kubota et al. 2001). 

LMXBs exhibit hard X-ray states which are very similar to the Black Hole Candidates ones and investigating the 
possible disc-jet coupling is very interesting and still matter of discussion. The power-law component could be produced 
by Componization by synchrotron emission in the jet (Bosch-Ramon et al. 2005, Fender 2004) as proposed for others 
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LMXBs: the power law component in Cir X-1 is explained by laria et al. (2002) as jet emission that has been resolved 
using radio interferometry by Fender & Kuulkers (2001); a significant correlation between radio and X-ray flux has 
been detected for 4U 1728-34 (Mighari et al. 2003), indicating a clear signature of disc-jet coupling. In our source, 
this hypothesis is strengthened by radio detection: Sydney University Molonglo Sky Survey catalog gives a flux of 7.5 
mJy at 843 GHz (Mauch et al. 2003). Assuming a radio spectral index of -0.5 to estimate the flux density at 5 GHz 
based on other observations, we found a radio loudness Pr/Px (as defined in Fender & Kuulkers 2001) in the range 
1.1-2.2 Jy/Crab, where Pr is the radio flux density at 5 GHz and Px is the peak X-ray flux measured in the soft band 
(< 12 keV), even though the two data set are not simultaneous. The radio loudness value is consistent with the one for 
Cir X-1 (Fender & Kuulkers 2001) and it is very high respect to other LMXBs and quite similar to that of Back Hole 
Candidates. These founding are in favour of the jet hypothesis as origin of the power-law observed for 4U 1636-53 
similarly to GX 5-1, Cir X-1 and GX 17+2 being all of them associated with variable radio sources (Fender & Hendry 
2000, and reference therein). 

Simultaneous high energy and radio observations during spectral transition are crucial to disentangle the power law 
emission in the LMXBs containing a neutron star. 

Finally, Laurent & Titarchuk (1999) suggest the detection of a power-law component at high energy to be a signature 
of presence black hole in an X-ray binary system. Our data are supporting result on GX 17+2 by Di Salvo et al. 
(2000), clearly showing this criterion proposed to distinguish black hole versus neutron star binaries is inadequate. 

We acknowledge the AST financial/programmatic support via contracts ASI-IR 046/04. A special thank to M. 
Federici and G. De Cesare for supervising the INTEGRAL data archive and software respectively. A very particular 
thank to K. Kretschmer for making data available before becoming public. 
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